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Abstract 

We recently proposed that, due to the top-quark-mass enhanced Yukawa coupling, the 
s-channel production of a charged scalar or pseudo-scalar from heavy quark fusion can 
be an important new mechanism for discovering non-standard spin-0 particles. In this 
work, we present the complete 0{a s ) QCD corrections to this s-channel production 
process at hadron colliders, and also the results of QCD resummation over multiple 
soft-gluon emission. The systematic QCD-improved production and decay rates at the 
FermiLab Tevatron and the CERN LHC are given for the charged top-pions in the 
topcolor models, and for the charged Higgs bosons in the generic two Higgs doublet 
model. The direct extension to the production of the neutral (pseudo-)scalars via bb 
fusion is studied in the minimal super symmetric standard model (MSSM) with large 
tan/3, and in the topcolor model with large bottom Yukawa coupling. 
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1. Introduction 

The top quark (t), among the three generations of fermions, is the only one with a large 
mass as high as the electroweak scale. This makes the top the most likely place to discover 
new physics beyond the Standard Model (SM). In a recent study two of us proposed that, 
due to the top-mass enhanced flavor mixing Yukawa coupling of the charm (c) and bottom 
(6) with a charged scalar or pseudo-scalar (0 ± ), the s-channel partonic process cb,cb — > ± , 
can be an important mechanism for the production of (f)^ at various colliders. From the 
leading order (LO) calculation M, we demonstrated that the FermiLab Tevatron Run-II has 
the potential to explore the mass range of the charged top-pions up to about 300-350 GeV 
in the topcolor (TopC) models || || . In this work, we compute the complete next-to-leading 
order (NLO) QCD corrections to the process qq' — > ± , which includes the one-loop virtual 
corrections and the contributions from the additional 0(a s ) processes, 

qq -> ± 5- and qg -> g> ± . (1) 

The decay width and branching ratio (BR) of such a (pseudo-) scalar are also included up to 
NLO to estimate the event rates. The QCD resummation of multiple soft-gluon radiation is 
also carried out, which provides a better prediction of the transverse momentum distribution 
of the (pseudo-)scalar particle. We shall choose the TopC model M as a benchmark of 
our analysis. The generalization to the generic type-Ill two-Higgs doublet model (2HDM) 
[fj. [| is straightforward since the Q CD-corrections are universal.[] The direct extension to the 
production of neutral (pseudo- )scalars via bb fusion is studied in the Minimal Supersymmetric 
SM (MSSM) || [7j with large tan/3 and in the TopC models with [/(l)-tilted large bottom 
Yukawa coupling [§, j3|. 

2. Charged Scalar Production via Charm-Bottom Fusion 

2.1. Fixed-Order Analysis up to 0(a s ) 

We study charged (pseudo-) scalar production via the top-mass-enhanced flavor mixing 
vertex c-b-^ The corresponding Yukawa coupling can be generally defined as ClL+CrR 
in which L = (1 — 75) /2 and R = (l + 7 5 )/2 . The total cross sections for the 4> + production 

1 We note that the finite part of the counter term to the q-q'-cjp' Yukawa vertex is renormalization-scheme- 
and model-dependent. 
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at hadron colliders (cf. Fig 1) can be generally expressed as 
a [h x h 2 -> <P + X) = f dasi T dx 2 [L/h^Q^fp/h^Q 2 ) + (a <-> (3)] o a(i {a$ - 

(2) 

where r = rn^/S, £1,2 = V^ e ±y , is the mass of \fS is the center-of-mass energy of 
the hifi2 collider, and f a /h{x,Q 2 ) is the parton distribution function (PDF) of a parton a 
with the factorization scale Q. The quantity a a ^ is the partonic cross section and has the 
following LO contribution for cb — > <p + (cf. Fig. la) [|TJ: 

where r = raff's with s~ the center-of-mass energy of the sub-process, and the terms sup- 
pressed by the small mass ratio (m^/m^) 2 have been ignored. Since we are interested in the 
inclusive production of the scalar 0, it is natural to choose the factorization scale Q to be its 
mass m<£, which is of O(10 2-3 ) GeV and much larger than the mass of charm or bottom quark. 
Hence, in this work, we will treat c and b as massless partons inside proton or antiproton 
and perform a NLO QCD calculation with consistent sets of PDFs || |[| |10| . 

The NLO contributions are of 0(a s ), which contain three parts: (i) the one-loop Yukawa 
vertex and quark self-energy corrections (cf. Fig. lb-d); (ii) the real gluon emission in the 
gff-annihilations (cf. Fig. le); (iii) s- and t-channel gluon-quark fusions (cf. Fig. lf-g). The 
Feynman diagrams coming from permutations are not shown in Fig. 1. Unlike the usual Drell- 
Yan type of processes (where the sum of the one-loop quark-wavef unction renormalization and 
vertex correction gives the ultraviolet finite result), we need to include the renormalization 
for the Yukawa coupling (yj) which usually relates to the relevant quark mass (Tn qj ), i.e., 
we have to add the counter term at the NLO (cf. Fig. Id) besides the contribution from 
the usual wavefunction renormalization Z qiq2 ^ = h(Z qi + Z q2 ) (cf. Fig.lc). This applies 
to the Yukawa interactions of the SM and MSSM Higgs bosons as well as the top-pions in 
the TopC models. It is clear that, for flavor-mixing vertex c-6-0 ± in the TopC model [cf. 
Eq. (|1I]) below], the counter-term of the Yukawa coupling is equal to the top quark mass 
counter-term Sm t /m t , which we determine from the top-quark mass renormalization in the 
on-shell scheme so that m t is the pole mass of the top quark. In other cases such as in the 
general 2HDM (type-Ill) || and the TopC models (with 6-Higgs or 6-pions) ||11|| , some of 
their Yukawa couplings are not related to quark masses or not of the above simple one-to- 
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one correspondence, and thus have their independent counter terms (Syj/yj). In addition 
to the virtual QCD-loop corrections, the contributions of the real gluon emission from the 
initial state quarks have to be included (cf. Fig. le). The soft and collinear singularities 
appeared in these diagrams are regularized by the dimensional regularization prescription at 
D = 4 — 2e dimensions. After summing up the contributions of virtual gluon-loop and real 
gluon-radiation (cf. Fig. lb-e), the ultraviolet and soft singularities separately cancel. The 
collinear singularities are still left over and should be absorbed into the renormalization of the 
PDF [fT2| . (The MS renormalization scheme is used in our calculation.) Finally, the gluon- 
quark fusion sub-processes (cf. Fig. lf-g) should also be taken into account and computed at 
general dimension-/}. All these results are separately summarized into the Appendix. 




(a) (b) (c) (d ) 




Figure 1: Representative diagrams for charged or neutral (pseudo-)scalar (dashed line) production 
from quark-antiquark and quark-gluon collisons at 0(oP s ) and 0(a).): (a) leading order contribution; 
(b-d) self-energy and vertex corrections (with counter term); (e) real gluon radiation in gg'-fusion; 
(f-g) s- and f-channel gluon-quark fusions. 

The hadron cross sections become regular after renormalizing the Yukawa coupling and 
the PDFs in (f2|), which are functions of the renormalization scale /i and the factorization scale 
Rf( = y/Q*)- The partonic NLO cross section ct^lo ( a ~^ + ^O contains the contributions 
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where r = ruff's and CV = 4/3. The mass counter term for the Yukawa vertex renormaliza- 
tion is determined in the on-shell scheme, i.e., 



5m t 
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-7 E + ln47r + ft 



(5) 



in the TopC model. Here, the bare mass m t o and the renormalized mass m t are related by 
m t0 = m t + 5m t and m t ~ 175 GeV is taken to be the top-quark pole mass. The finite part 
of the mass counter term is ft = 3 In \ji 2 /mf\ +4 in the TopC model, where ft > for 
/i > mje _2//3 ~ 90 GeV. In the following, we shall choose the QCD factorization scale [if 
(set as the invariant mass yQ^) and the renormalization scale /i to be the same as the scalar 
mass, i.e., a/Q 2 = /i = m^,, which means that in (|j) the factor In (m 2 /<5 2 ) vanishes and the 
quantity ft becomes 

ft = 3 In [my ml] + 4 . (6) 

For the case of m^ 3> m t , the logarithmic term In (m 2 /m 2 ) becomes larger for m^ ^> m t , 
and its contributions to all orders in a s In (m^/m 2 ) may be resummed by introducing the 
running Yukawa coupling t/t{f^), or correspondingly, the running mass m t (fi). In the above 
formula, m t is the pole mass (mf o1 ~ 175 GeV) and is related to the one- loop running mass 
via the relation Jl3| 



m t (fi) = m t (m 



r 1 ) 



3C F , s, ^ 
1 : — a s (fi) In 



2 i 
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(7) 



Using the renormalization group equation, one can resum the leading logarithms to all orders 
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TM and obtains 



m t (fi) = mt(m 



a s (n) 



9C F 



(8) 



a s {m p t oX )_ 

with rif = 6 for li > m t . Thus, to include the running effect of the Yukawa coupling, we 
can replace the (mf ol ) 2 -f actor (from the Yukawa coupling) inside the square of the S'-matrix 
element [up to 0(a s )} by the running factor 

2" 
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(9) 



where the logarithmic term in the bracket [• • ■] is added to avoid double-counting with the 
resummed logarithms inside m 2 (fi). It is clear that this [1 + (Cpa s (li) /2tt) O] factor will 
cancel the f2-term inside the NLO hard cross section Aa cb - in Eq. (^) at 0(a s ), so that the net 
effect of the Yukawa vertex renormalization (after the resummation of leading logarithms) is 
to replace the relevant tree-level on-shell quark mass (related to the Yukawa coupling) by its 
MS running mass [cf. Eq. (|^)] and remove the Jl-term in Eq. (||). When the physical scale 
fj, (chosen as the scalar mass m^) is not much larger than m t , the above running effect is 
small since the ln(/i/m t ) factor in the Yukawa counter-term 5mt/mt is small. However, the 
case for the neutral scalar production via the bb annihilation can be different. When the loop 
correction to the (fP-b-b Yukawa coupling contains the logarithm ln(/i/m&), which is much 
larger than ln(/i/m 4 ), these large logarithms should be resummed into the running coupling, 
as we will do in Section 4. 

In the TopC model, there are three pseudo-scalars, called top-pions, which are predicted 
to be light, with a mass of 0(100 ~ 300) GeV. The relevant Yukawa interactions for top-pions, 
including the large £_r-cr flavor-mixing, can be written asQ Jl] 

m t tan/3 r tt tt *_ ./K-r^tt* j^bb i— h + . 

iK UR K UL t L t R n t + V2K UR K DL t R b L TT t + 

v L (10) 

iK^K^IEcn^ + y/2K%*K» L cRb L «t + h.c] , 

where tan/3 = \J (v/vt) 2 — 1 ~ 0(4—1.3) with the top-pion decay constant fi~O(60— 150) GeV, 
and the full vacuum expectation value (VEV) v ~ 246 GeV (determined by the Fermi con- 
stant). The analysis from top-quark decay in the Tevatron ti events sets a direct lower bound 

2 As pointed out in Ref. an important feature deduced from (10) is that the charged top-pion irf 1 mainly 
couples to the right-handed top (tn) or charm (cr) but not the left-handed top (Jul) or charm (cl), in contrast 
to the standard W-t-b coupling which involves only th- This makes the top-polarization measurement very 
useful for further discriminating the signal from the background events. 
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on the charged top-pion mass to be larger than about 150 GeV 1TR [21. The existing low en- 



ergy LEP/SLD measurement of Rb, which slightly lies above the SM value by about 0.9a [16 



also provides an indirect constraint on the top-pion Yukawa coupling = (\/2m t /v)taa/3 
due to the one-loop contribution of charged top-pions to Rb- However, given the crude ap- 
proximation in estimating the top-pion loops (with all higher order terms ignored) and the 
existence of many other sources of contributions associated with the strong dynamics, the 
indirect Rf, constraint is not conclusive 0. For instance, it was shown that the 3a Rb bound 
from the one-loop top-pion correction can be fully removed if the top-pion decay constant 
v t is increased by about a factor of 2 (which is the typical uncertainty of the Pagels-Stokar 
estimate) |2|, |T7[|; also, the non-perturbative contributions of the coloron-exchanges can shift 
the Rb above its SM value and tend to cancel the negative top-pion corrections. Due to 
these reasons, it is clear that the inconlusive ii^-bound in the TopC models should not be 
taken too seriously. Nevertheless, to be on the safe side, we will impose the roughly estimated 
-Rb-constraint in our current analysis of the TopC model, by including only the (negative) 
one-loop top-pion contribution as in Ref. |17|].P| As shown in Fig. [2]a, the current 3a i?f,-bound 
requires a smaller top-pion Yukawa coupling, C R ~ 1.3 — 2 (or, tan/3 ~ 1.3 — 2), for the low 
mass region of ~ 200 — 500 GeV. Since the top-pion decay constant Vt is related to tan/3, 
this also requires v t to be around 150 ~ 100 GeV for m n ± ~200 — 500 GeV (cf. Fig. 0b). For 
comparison, the usual Pagels-Stokar estimate of Vt (by keeping only the leading logarithm 
but not constant terms), v\ = (N c /87r 2 )m 2 hiA 2 /m 2 , gives t> < ~64 — 97 GeV for the topcolor 
breaking scale A ~ 1 — 10 TeV, where a typical factor of 2 ~ 3 uncertainty in the calculation 



of v 2 is expected H, 18]. This estimate is slightly lower than the -R&-constrained values of Vt 



in Fig. [2b, but is still in reasonable consistency (given the typical factor of v2~\/3 error in 

the leading logarithmic Pagels-Stokar estimate of Vt). 

In (|T0"D, Kjjl,r and Kbl,r are defined from diagonalizing the up- and down-type quark 

mass matrices M v and M D : K ] UL MuK UR = M$ ia , K ] DL M D K DR = M^ ia , with M$ ia = 

diag(m u , m c , m t ) and M^ ia = diag(m c ;, m s , m^). For the class-I TopC models [II], we have 

constructed || a realistic and attractive pattern of Km and Kol so that the well-constrained 

3 However, it is important to keep in mind that such a rough i?6-bound is likely to over-constrain the 
top-pion Yukawa coupling since only the negative one-loop top-pion correction (but nothing else) is included 
in this estimate. A weaker i?{,-bound will less reduce the top-pion Yukawa coupling and thus allow larger 
production rates of charged top-pions at colliders which can be obtained from our current analysis by simple 
re-scaling. 
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Cabibbo-Kobayashi-Maskawa (CKM) matrix V (= K ul Kdl) can be reproduced in the 
Wolfenstein-parametrization |19| and all potentially large contributions to the low energy 
data (such as the K-K, D-D and B-B mixings and the b — > 37 rate) can be avoided U . We 
then found that the right-handed rotation matrix Kjj R is constrained such that its 33 and 32 
elements take the values as M 



0.99-0.94, K* R <y/l-K$ R '~ 0.11-0.33, (11) 

which show that the t R -c R flavor mixing can be naturally around 10 — 30%. 

For the current numerical analysis we consider a benchmark choice |l[ based upon the 
above TopC model: 

Cf = C^{R b constrained), = CrK% r ~ 0% x 0.2, Cf = Cf = 0. (12) 

It is trivial to scale the numerical results presented in this paper to any other values of 



Cl,r when needed. Unless specified otherwise, we use CTEQ4M PDF [2D| to calculate the 
rates. Note that CTEQ4M PDFs are consistent with the scheme used in the current study 
which treats the initial state quarks as massless partons in computing the Wilson coefficient 
functions. The only effect of the heavy quark mass is to determine at which scale Q this 
heavy quark parton becomes active. f\ In our case, the scale Q = ^> m c , mi,. 

In Fig. H|, we present the total cross sections for the charged top-pion production as 
functions of its mass, at the Tevatron (a pp collider at 1.8 and 2TeV) and the LHC (a pp 
collider at 14TeV). We compare the improvements by including the complete NLO results 
[cf. (|])] and by including the resummed running Yukawa coupling or running mass [cf. (|8|)]. 
For this purpose, we first plot the LO total cross sections with the tree-level Yukawa coupling 
[dash-dotted curves, cf. (^|) and (|T2|)1 and with the resummed running Yukawa coupling or 
running mass [dotted curves, cf. (|3]) and (H)]; then we plot the NLO cross sections with the 
one-loop Yukawa coupling [dashed curves, cf. (|j)] and with the resummed running Yukawa 
coupling or running mass [solid curves, cf. (£|), (§) and ([])] . We see that at the LHC there is 
a visible difference between the pure LO results with tree-level Yukawa coupling (dash-dotted 
curves) and other NLO and/or running- coupling improved results. But at the Tevatron, the 
LO results with running Yukawa coupling (dotted curves) are visibly smaller than the results 
4 This is the Collins- Wilczek-Zce (CWZ) scheme §1. 
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in all other cases for > 300 GeV. This shows that without the complete NLO calculation, 
including only the running Yukawa coupling in a LO result may not always warrant a better 
improvement. Finally, the comparison in Fig. [3] shows that the resummed running Yukawa 
coupling or top-mass [cf. Eq. @] does not generate any significant improvement from the 
one-loop running. This is because the top-mass is large and a s In (m^/mf) is small for up 
to 1 TeV. Thus, the improvement of the resummation in (§) has to come from higher order 
effects of a s ln (m^/m^). However, as to be shown in Sec. 4, the situation for summing over 
powers of a s In (m 2 ,/m%) is different due to m b -C m t , . 

Fig. f| is to examine the individual NLO contributions to the charged top-pion production 
via the qq' and qg sub-processes, in comparison with the full NLO contributions.^ The LO 
contributions are also shown as a reference.^ [Here q denotes the heavy charm or bottom 
quark.] In this figure, there are three sets of curves for the charged top-pion production cross 
sections: the highest set is for the LHC (y/S = 14 TeV), the middle set is for the upgraded 
Tevatron (y/S = 2 TeV), and the lowest set is for the Tevatron Run I (y/S = 1.8 TeV). The 
LO cross sections are plotted as dotted lines while the NLO cross sections as solid ones. 
The dashed lines show the contributions from the gg'-fusion sub-processes, and the dash- 
dotted lines describe the contributions from the gg-fusion sub-processes. The gg-fusion cross 
sections are negative and are plotted by multiplying a factor of —1, for convenience. For 
a quantitative comparison of the individual NLO contributions versus the full NLO results, 
we further plot, in Fig. [|, the ratios (called ^-factors) of the different NLO contributions 
to the LO cross section by using the same set of CTEQ4M PDFs. The solid lines of Fig. |5| 
show that the overall NLO corrections to the pp,pp — > ± V processes are positive for 
above ~150 (200) GeV and lie below ~15 (10)% for the Tevatron (LHC) in the relevant mass 
region. This is in contrast with the NLO corrections to the boson production at hadron 
colliders, which are always positive and as large as about 25% at the Tevatron |[22|| . The 

T qq , and &CTqg,gq 
Tqq> 

the case of ± production the size of negative Aa qg ^q piece becomes comparable with that 

of the positive Aa q q> such that a non-trivial cancellation occurs. 

Unless specified, qg includes both qg and qg contributions. 
6 With the exception of Figs. ||, || and [l^, we only show our numerical results with the resummed 
running Yukawa coupling or running mass. 



reason of this difference originates from the differences in the Aa qq -> and Aa qg ^q for (p^ and 
production. While in the case of production the positive Aa q q> piece dominates, in 
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While it is reasonable to take the renormalization and the factorization scales to be for 

predicting the inclusive production rate of <ft + , it is desirable to estimate the uncertainty in 

the rates due to different choices of PDFs. For that purpose, we examine a few typical sets of 

PDFs from CTEQ4, which predict different shapes of charm, bottom and gluon distributions. 

As shown in Table |l] and Fig. |6|, the uncertainties due to the choice of PDF set are generally 

within ±20% for the relevant scalar mass ranges at both the Tevatron and the LHC. 

Table 1: Cross sections in fb for charged top-pion production in the TopC model at the upgraded 
Tevatron and the LHC are shown, by using four different CTEQ4 PDFs. They are separately given 
for the LO and NLO processes, and for the qq — > (p + X and qg — > (j) + X sub-processes. At the 
upgraded Tevatron the top number is for = 200 GeV, the middle is for = 300 GeV, and the 
bottom is for = 400 GeV. At the LHC the top number is for = 400 GeV, the middle is for 
= 700 GeV, and the lowest is for = 1 TeV. 



Collider 


Upgraded Tevatron (2 TeV) 


LHC (14 TeV) 


Process \ PDF 


4A1 


4M 


4A5 


4HJ 


4A1 


4M 


4A5 


4HJ 




367 


382 


376 


387 


5380 


5800 


6060 


5890 


LO 


42.6 


43.7 


41.5 


46.6 


863 


901 


896 


906 




6.88 


7.05 


6.56 


8.38 


235 


240 


232 


241 




370 


402 


412 


407 


5430 


6080 


6510 


6170 


NLO 


45.6 


48.6 


47.9 


51.6 


912 


976 


997 


981 




7.70 


8.21 


7.89 


9.56 


255 


266 


264 


268 




551 


584 


585 


590 


7530 


8290 


8740 


8400 


qq — > (j) + X 


64.5 


67.4 


65.5 


71.7 


1210 


1280 


1290 


1290 




10.6 


11.1 


10.5 


13.0 


331 


341 


335 


343 




- 180 


- 181 


- 174 


- 183 


-2100 


-2200 


-2240 


-2230 


qg - 0+X 


-19.2 


-18.9 


-17.5 


-19.9 


-299 


-302 


-293 


-303 




-2.94 


-2.86 


-2.59 


-3.34 


-76.0 


-74.7 


-70.6 


-75.0 



2.2. Analysis of Multiple Soft-Gluon Resummation 

The a s corrections to the (pseudo-)scalar production involve the contributions from the 
emission of virtual and real gluons, as shown in Figs. 1(b), (c) and (e). As the result of 
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the real gluon radiation, the (pseudo-) scalar particle will acquire a non- vanishing transverse 
momentum (Qt)- When the emitted gluons are soft, they generate large logarithmic con- 
tributions of the form (in the lowest order): a s ln m (Q 2 /Q T ) /Qti where Q is the invariant 
mass of the (pseudo-)scalar, and m — 0, 1. These large logarithms spoil the convergence of 
the perturbative series, and falsify the 0(a s ) prediction of the transverse momentum when 
Qt < Q- 

To predict the transverse momentum distribution of the produced (pseudo-) scalar, we uti- 
lize the Collins-Soper-Sterman (CSS) formalism P5| , resumming the logarithms of the type 
a™ ln m {Q 2 /Q T ) jQ\ , to all orders n in a s (m = 0, 2n — 1). The resummation calculation 



is performed along the same line as for vector boson production in Ref. ||22|| . Here we only 
give the differences from that given in Ref. [P^j . But for convenience, we also list the A^ l \ 
A^ 2 \ and coefficients of the Sudakov exponent, which have been used in the current 
analysis: 

AW(d) = C F , B<»(C 1 = b ,C 2 = l) - 



'2 Cf 



A™ (Ci 



67 7T 2 \ 



(13) 



v 36 12 ) N ° 18 Hf _ ' 
where Cp = 4/3 is the Casimir of the fundamental representation of SU(3), Nq = 3 is the 
number of SU(3) colors, and n/ is the number of light quark flavors with masses less than 
Q. In the above we used the canonical values of the renormalization constants C\ = bo, and 
C 2 = 1. 

To recover the O (a s ) total cross section, we also include the Wilson coefficients Cj£\ 
among which C^j differs from the vector boson production (here i denotes quark or antiquark 
flavors, and a = qi or gluon g). Explicitly, 

M,/i,CVC 2 ) = 5 jk 8(l-z), C^iz^^Ci/d) 

c®{z&frC x /c*) 



0. 



cV(z,b,n,Ci/C 2 ) 



+ 6(1 -z) 



C h 



In 



fib 



In 2 



ln[^ 



b C 2 



3/4 



V 



+ 



4 



9 
16 



bo 



(14) 



where P„-./L is the 0(a s ) gluon splitting kernels [24, given in the Appendix. In the above 
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expressions, V = Vdy = — 8 + 7r 2 for the vector boson production p2"fl , and V = V$ = n 2 for 
the (pseudo-) scalar production, when using the running mass given in Eq. (§) for the Yukawa 
coupling. Using the canonical values of the renormalization constants, ln(/i6/6 ) vanishes, 
because /i = C\/b = b /b. 

The only remaining difference between the resummed formulae of the vector boson and 
(pseudo-) scalar production is in the regular (Y) terms, which comes from the difference of 
the O (a s ) real emission amplitude squares (cf., the definitions of Tz 1 and T^ 1 in Appendix 
C of Ref. |22|] and Eqs. fl22|) and ( p5|) of this paper). The non-perturbative sector of the CSS 
resummation (the non-perturbative function and the related parameters) is assumed to be 
the same as that in Ref. |22|| . 

As described in Ref. [p2 |, the resummed total rate is the same as the 0(a s ) rate, when 
we include and and switch from the resummed distribution to the fixed order one 
at Qt = Q- When calculating the total rate, we have applied this matching prescription. In 
the case of the (pseudo-) scalar production, the matching takes place at high Q T ~ Q values, 
and the above matching prescription is irrelevant when calculating the total rate because the 
cross sections there are negligible. Thus, as expected, the resummed total rate differs from 
the 0(a s ) rate only by a few percent. Since the difference of the resummed and fixed order 
rate indicates the size of the higher order corrections, we conclude that for inclusive (pseudo- 
scalar production the 0(a 2 ) corrections are likely much smaller than the uncertainty from 
the parton distribution functions (cf. Fig.^j) . 

In Fig. [7|, we present the numerical results for the transverse momentum distributions of 
the charged top-pions (in TopC model) and the charged Higgs bosons (in 2HDM) produced 
at the upgraded Tevatron and the LHC. The solid curves show the resummation prediction 
for the typical values of m^. The dashed curves, from the 0(a s ) prediction, are irregular as 
Qt — » 0. The large difference of the transverse momentum distributions between the results 
from the resummation and fixed-order analyses throughout a wide range of Qt shows the 
importance of using the resummation prediction when extracting the top-pion and Higgs 
boson signals. We also note that the average value of Qt varies slowly as increases and 
it ranges from 35 to 51 GeV for between 250 and 550 GeV at the 14 TeV LHC, and from 
23 to 45 GeV for between 200 and 300 GeV at the 2 TeV Tevatron. 
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Figure 2: Estimated current 3<r-bounds in the TopC model and 2HDM-III: (a) the 3a upper bound 
on the top-pion Yukawa coupling C^; (b) the 3a lower bound on the top-pion decay constant; 
(Here, in (a) and (b), the solid curves are derived from the combined LEP/SLD data of R^ xp = 
0.21656 ± 0.00074 while dashed curves are from the same 3a combined experimental error but with 
the central -Revalue equal to Rf M = 0.2158.) (c) the ^-predictions of 2HDM-III with coupling 
^tt = 1-0 and 1.5 (solid curves) and the 3a ^-bounds (dashed lines). 



12 



Top-pion Production 
(CTEQ4M) 




1 00 



200 300 400 500 



600 



700 



800 900 1000 

m (GeV) 



Figure 3: Top-pion production cross sections at the present Tevatron, upgraded Tevatron, and 
the LHC. For each collider we show the NLO cross section with the resummed running Yukawa 
coupling (solid), and with one-loop Yukawa coupling (dashed), as well as the LO cross section with 
resummed running Yukawa coupling (dotted) and with tree-level (dash-dotted) Yukawa coupling. 



3. Hadronic Decays of Charged Scalars to 0(a s ) 

In the TopC models, the current Tevatron data from the top quark decay into charged 
top-pion (irf) and 6-quark already requires the mass of irf 1 to be above ~150 GeV ||, [15| . In 
the current analysis, we shall consider m 7Tt > m t + rrib, so that its dominant decay channels 
are — > tb,cb . 

The decay width of nf (= 0*), including the 0(a s ) QCD corrections, is given by [EB, j27|: 
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in which Q = a/Q 2 is the invariant mass of (f) ± . The small bottom and charm masses are 
ignored so that r = (m t /m^) 2 for tb final state and r = for c6 final state. Thus, for 
^ — ► c6, the quantity 7£ reduces to 7?. = 17/2 — VL . In Fig. || we present the results for total 
decay widths of <p + and branching ratios of <p + — > t6 in the TopC model and the 2HDM. For 
the 2HDM, we also show the branching ratios of the W + h° channel, which is complementary 
to the tb channel. The NLO (solid) and LO (dashed) curves differ only by a small amount. In 
the same figure, the i^-factor, defined as the ratio of the NLO to the LO partial decay widths, 
is plotted for the + — > tb (solid) and — > cb (dashed) channels. Here, the sample results for 
the 2HDM are derived for the parameter choice: a = and (M h , M A ) = (100, 1200) GeV. 

With the decay width given above, we can study the invariant mass distribution of tb for 
the s-channel + -production: 

(Q 2 r /m^) Br [0+^ tb] 



da 
~dQ 2 



[h x h 2 



-X)^tbX] = a [h l h 2 -^(j) + {Q)X] 



71 



{Q 2 -miY + (Q^^y 



(16) 



tb. 



where and Br [</> + ^ tb\ are the total decay width of <p + and the branching ratio of (p + 
respectively, which are calculated up to the NLO. We note that the one- loop box diagrams 
with a virtual gluon connecting the initial state quark and final state quark (from the hadronic 
decay of <p) have vanishing contribution at 0(a s ) because the scalar <p is color- neutral. In 
Fig. |9|a and Fig. |i~0|a, we plot the invariant mass distribution for t-b and t-b pairs from (p^ 
(top-pion signal) and iy ± * (background) decays in the TopC model. In these plots, we 
have included the NLO contributions, as a function of Q, to the W^* background rate at 
the Tevatron and the LHC. The overall i^-factor (after averaging over the invariant mass 
Q) including both the initial and final state radiations is about 1.4 (1.34) for the Tevatron 
(LHC) @. The total rate of up to the NLO is about 0.70 [0.86] pb and 11.0 pb at the 
1.8 [2] TeV Tevatron and the 14 TeV LHC, respectively. 
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Figure 4: Cross sections for the charged top-pion production in the TopC model at the present 
Tevatron, upgraded Tevatron and the LHC. The NLO (solid), the qq' (dashed) and qg (dash-dotted) 
sub-contributions, and the LO (dotted) contributions are shown. Since the qg cross sections are 
negative, they are multiplied by —1 in the plot. The cross sections at y/S = 1.8 TeV are multiplied 
by 0.1 to avoid overlap with the y/S = 2 TeV curves. 
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Figure 5: The K- factors for the 4> + production in the TopC model are shown for the NLO (K 
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Figure 6: The ratios of NLO cross sections computed by four different sets of CTEQ4 PDFs relative 
to that by the CTEQ4M for charged top-pion production at the upgraded Tevatron (a) and the 
LHC (b). 
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Figure 8: Total decay widths of (p+ and BRs of cf> + -> tb in the TopC model and 2HDM. (For the 
2HDM, the BR of the W + hP channel is also shown, which is complementary to the tb channel.) In 
Fig. (a) and (b), the NLO (solid) and LO (dashed) curves differ only by a small amount. In Fig. 
(c), the if-factor, which is defined as the ratio of the NLO to the LO partial decay widths, is shown 
for the 4> + — > tb (solid) and — > cb (dashed) channels. The sample results for the 2HDM in this figure 
are derived for the parameter choice f^ c ) = (1.5, 1.5), a = 0, and (m^, uia) = (120, 1200) GeV. 



19 




180 200 220 240 260 280 300 320 340 360 



Q (GeV) 

Figure 9: Invariant mass distribution of t-b and i-b pairs from ^ (signal) and W^* (background) 
decays at the Tevatron Run-II for the TopC model (a), and 2HDM with Higgs mixing angles a = 
(b), and a = tt/2 (c). We show the signal for = 200, 250, 300 and 350 GeV. The solid curves 
show the results from the NLO calculation, and the dashed ones from the LO analysis. 
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Figure 10: Invariant mass distribution of t-b and i-b pairs from ± (signal) and W^* (background) 
decays at the LHC for the TopC model (a), and for the 2HDM with the Higgs mixing angles a = 
in (b), and a = tt/2 in (c). Here the charged pseudo-scalar or scalar mass are chosen as the typical 
values of = 250, 400 and 550 GeV. The solid curves show the results by the NLO calculation, 
while the dashed ones come from the LO analysis. 
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Figure 11: Invariant mass distributions of W + -h° and W~-h° pairs from cjy^ (s-channel resonance) 
and W^* (s-channel non-resonance) decays at the Tevatron Run-II, and at the LHC, for the 2HDM 
with Higgs mixing angles a = 0. We show the signal for = 200, 250 and 300 GeV at the Tevatron 
(a), and for = 250, 400 and 550 GeV at the LHC (b). The solid curves show the results of the 
NLO calculation, and the dashed ones of the LO analysis. 



Before concluding this section, we discuss how to generalize the above results to the generic 
2HDM (called type-Ill in which the two Higgs doublets $i and $2 couple to both up- 
and down-type quarks and the ad hoc discrete symmetry [2~9| is not imposed. The flavor- 
mixing Yukawa couplings in this model can be conveniently formulated under a proper basis 
of Higgs doublets so that <$i>= (0,v/V~2) T and <<J?2>= (0,0) T . Thus, the diagonalization 
of the fermion mass matrix also diagonalizes the Yukawa couplings of $1, and all the flavor- 
mixing effects are generated by Yukawa couplings (YV and Y?) of $2 which exhibit a natural 
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hierarchy under the ansatz j|, [5| 



Y; 



U,D 



(17) 



with £^- D ~ 0(1). This ansatz highly suppresses the flavor-mixings among light quarks 
and identifies the largest mixing coupling as the one from the t-c or c-t transition. A recent 
renormalization group analysis |3(J shows that such a suppression persists at the high energy 
scales. The relevant Yukawa interactions involving the charged Higgs bosons H ± are ffl: 



t R [%V) th b L - t L (VY D ) tb b R +cK (Y^V) C b b L -cE (VY D ) cb b R 



+ h.c. 



t R Y t ^b L + cRY tc 



u* 



+ h.c. + (small terms) 



where Y# = £^ x (y/2m t /i 



and Y? 



tc 



(18) 

£g x (^2m t m c /v ) ~ x 9% , in which 
£,Y C ~ 0(1) is allowed by the current low energy data |52]| . As a result, the Yukawa counter 
term in Fig. Id involves both 5m t and 5m c . Consequently, we need to replace the NLO 
quantity Q in the finite part of the Yukawa counter term [cf. the definition below (|])] by 



fi(2HDM) = 3 In [ml/(m t m c )] +4 



(19) 



for the type-Ill 2HDM. In the relevant <f> -c-b coupling of this 2HDM, we note that, similar 
to the case of the TopC model, only the right-handed charm is involved i.e., 



nth 



Ccb 
L 



0. 



ct * eg x 9%. 



(20) 



where the parameters (£ t ^, £g) are expected to be naturally around 0(1). We have examined 

the possible constraint of £^ from the current R b data and found that the values of £^ ~ 1.0— 1.5 

are allowed for ra#± > 200 GeV (cf. Fig.|2|c).[] The production cross section of H ± in this 

2HDM can be obtained by rescaling the result of the TopC model according to the ratio of 

the coupling-square [C^(2HDM)/C£r(TopC)] 2 ~ [0.09f£/0.2C#(TopC)] 2 (which is about 1/7 

for £Yc = 1-5 and the charged scalar mass around 400 GeV). 

Finally, we note that there are three neutral Higgs bosons in the 2HDM, the CP-even 

scalars (h°, H°) and the CP-odd pseudo-scalar A . The mass diagonalization for h° and H° 

7 Our calculation of i?& in the 2HDM-III is consistent with those in Ref. |3l| and Ref. |52| after using the 
same inputs. Note that a larger value of £^t than ours was chosen for the solid curve in Fig. 3 of Ref. p2| . 
We thank L. Reina for clarifying the inputs of Ref. and for useful discussions. 
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induces the Higgs mixing angle a. The low energy constraints on this model require |5], |32 
m hi rrtu — m H ± — m A or < mjj± < rnn- For the case of m#± > m h o + Myy, 
the if 1 * 1 — > W ± h° decay channel is also open. Taking, for example, a = and (m^, m^) = 
(120, 1200) GeV, we find from Fig. |8|b that the tb and Wh° decay modes are complementary 
at low and high mass regions of the charged Higgs boson H ^ . In Figs. |]b-c and |lO|b-c, we plot 
the invariant mass distributions of tb and tb pairs from if ± (signal) and W (background) 
decays in the 2HDM at the 2TeV Tevatron and the 14TeV LHC, with the typical choice 
of the parameters: (fg, f£) = (1.5, 1.5) in Eq. (gU|), (m h , m A ) = (120, 1200) GeV, and 
a = or 7r/2. [A larger value of will simutaneously increase (reduce) the BR of tb (Wh°) 
mode.] We see that, due to a smaller c-b-H ± coupling [cf. (|20|)1, it is hard to detect such a 
charged Higgs boson with mass m H ± > 250 GeV at the Tevatron Run-II. We then examine 
the potential of the LHC for the high mass range of H ± . Similar plots are shown in Figs. [5|b-c 
for a = and a = tt/2, respectively. When cos a is large (e.g., a = 0), the branching ratio of 
the tfo-channel decreases as m H ± increases (cf. Fig. |8|b), so that the LHC does not significantly 
improve the probe of the large mn±, range via the single-top mode (cf. Fig. |T0"|b). In this 
case, the W ± h° channel, however, becomes important for large mn±, as shown in Fig. |ll] 
(cf. Fig. [8]b, for its decay branching ratios) since the H^-W^-hP coupling is proportional to 
cos a p. On the other hand, for the parameter space with small cos a (e.g., a = tt/2), the 
W^hP channel is suppressed so that the single-top mode is important even for large mass 
region of H .| This is illustrated in Fig. |lO| c at the LHC for a = tt/2. In order to probe the 
whole parameter space and larger mn±, it is important to study both tb and Wh° (or WH°) 
channels. 

4. Generalization to Neutral Scalar Production via bb Fusion 

The QCD corrections are universal so that the generalization to the production of neutral 
scalar or pseudo-scalar 0° via the bb fusion is straightforward, i.e., we only need to replace 
©by 

Vl((f) bb) = 3 In [ml/m 2 b ] + 4 , (21) 

in which is the mass of 0°. The finite piece of the Yukawa renormalization [cf. the quantity 

8 Note that the H^-W^-H coupling is proportional to sin a and is thus enhanced for small cos a. In this 
case, the WH° mode may be important provided that H° is relatively light. We will not further elaborate 
this point here since it largely depends on the mass of H°. 
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f2 in (|5p] is scheme-dependent. We can always define the (jP-b-b Yukawa coupling as y/2mi,/v 
times an enhancement factor K so that the Yukawa counter term is generated by Smb/mb-f\ 
After resumming the leading logarithmic terms, [a s ln(m|/m^)] n , via the renormalization 
group technique, the net effect of the Yukawa renormalization is to change the Yukawa 
coupling or the related quark-mass into the corresponding MS running coupling or mass, as 
discussed in the previous section. 

The bb decay branching ratios of the neutral Higgs bosons in the MSSM with large tan/3 



are almost equal to one [34]. The same is true for the 6-Higgs or 6-pion in the TopC model 

@. It has been shown that at the Tevatron, the bb dijet final states can be properly identified 

p5| . The same technique developed for studying the resonance of the coloron or techni-p in 

the bb decay mode can also be applied to the search of the neutral Higgs bosons with 

large bottom Yukawa coupling. When the neutral scalar or pseudo-scalar </>° is relatively 

heavy, e.g., in the range of 0(250 — 1000) GeV, the QCD dijet backgrounds can be effectively 

removed by requiring the two 6-jets to be tagged with large transverse momenta (Pt) because 

the Pt of each 6-jet from the 0° decay is typically at the order of Hence, this process 

can provide complementary information to that obtained from studying the cjPbb associate 

production |36|, |33|, |37j . 

We first consider the production of the neutral Higgs boson 0°, which can be either 

A , h°, or if , in the MSSM with large tan/?, where the corresponding Yukawa couplings 

to bb and r + r~ are enhanced relative to that of the SM since yo/y^j 1 is equal to tan/3, 

— sin a/ cos/?, or cos a/ cos/3, respectively, at the tree-level. In the large tan/3 region, the 

MSSM neutral Higgs bosons dominantly decay into bb and r + r~ final states, which can be 

detected at the hadron colliders. In comparison with the recent studies on the <f)°bb and 

0°t + t~ associate production, we expect the inclusive 0° production via the 66-fusion 

would be more useful for being relatively heavy (e.g., > 200 — 300 GeV) because 

of the much larger phase space as well as a better suppression of the backgrounds in the 

high Pt region. The total LO and NLO cross sections for the inclusive production process 

pp,pp — > A°X at the Tevatron and the LHC are shown in Figs. |T2] a and b, in parallel to 

Figs. and § for the case of charged top-pion production. Here, we have chosen tan/3 = 40 

9 This specific definition works even if the Yukawa coupling is not related to any quark mass. For instance, 
the bottom Yukawa couplings of the 6-Higgs and 6-pion in the TopC model ^ |33| are independent of quark 
masses because the 6-Higgs does not develop VEV. 
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for illustration. The cross sections at other values of tan/3 can be obtained by multiplying the 



scaling factor (tan/3/40) . From Fig. |l2|a, we see a significant improvement from the pure 
LO results (dash-dotted curves) by resumming over the large logarithms of m|/m^ into the 
running Yukawa coupling. The good agreement between the LO results with running Yukawa 
coupling and the NLO results is due to a non-trivial, and pro cess- dependent, cancellation 
between the individual 0(a s ) contributions of the bb and bg sub-processes. In contrast to the 
production of the charged top-pion or Higgs boson via the initial state cb or cb partons, the 
neutral Higgs boson production involves the bb parton densities. The i^-factors for the ratios 
of the NLO versus LO cross sections of pp/pp — > A°X are presented in Fig. [13| for the MSSM 
with tan/3 = 40. The main difference is due to the fact that the individual contribution by 
the 0(a s ) 6^-fusion becomes more negative as compared to the case of the charged top-pion 
production shown in Fig. |5|. This makes the overall .ff -factor of the NLO versus LO cross 
sections range from about — (16~17)% to +5% at the Tevatron and the LHC. In parallel to 
Table 1 and Fig. [], we have examined the uncertainties of the CTEQ4 PDFs for the A°- 
production at the Tevatron and the LHC, and the results are summarized in Table [2] and 
Fig. |l^). We also note that, similar to the charged Higgs boson production, the resummed 
total rate for the neutral Higgs boson production is not very different from its NLO rate. 

The transverse momentum (Qt) distributions of A , produced at the upgraded Tevatron 
and at the LHC, are shown in Fig. [15] for various A masses (m^) with tan/? = 40. The 
solid curves are the result of the multiple soft-gluon resummation, and the dashed ones are 
from the 0(a s ) calculation. The shape of these transverse momentum distributions is similar 
to that of the charged top-pion (cf. Fig. [F|). The fixed order distributions are singular as 
Qt — > 0, while the resummed ones have a maximum at some finite Qt and vanish at Qt = 0. 
When Qt becomes large, of the order of ttla, the resummed curves merge into the fixed order 
ones. The average resummed Qt varies between 25 and 30 (40 and 60) GeV in the mass 
range of m A from 200 to 300 (250 to 550) GeV at the Tevatron (LHC). 

We also note that for large tan/?, the SUSY correction to the running <j)°-b-b Yukawa 
coupling is significant and can be included in a way similar to our recent analysis of the 
<f)°bb associate production ||33[| . To illustrate the SUSY correction to the b- Yukawa coupling, 
we choose all MSSM soft-breaking parameters as 500 GeV, and the Higgs mixing parameter 
fi = ±500 GeV. Depending on the sign of fx, the SUSY correction to the (fi°-b-b coupling can 
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either take the same sign as the QCD correction or have an opposite sign ||33|| . In Fig. |12|c, 
the solid curves represent the NLO cross sections with QCD correction alone, while the 
results including the SUSY corrections to the running bottom Yukawa coupling are shown 
for /i = +500 GeV (upper dashed curves) and fi = — 500 GeV (lower dashed curves). As 
shown, these partial SUSY corrections can change the cross sections by about a factor of 
2. The above results are for the inclusive production of the CP-odd Higgs boson A in the 
MSSM. Similar results can be easily obtained for the other neutral Higgs bosons (h° and 
H°) by properly rescaling the coupling strength. We also note that in the large tan/3 region, 
there is always a good mass-degeneracy between either h° and A (in the low mass region 
with rriA ~ 120 GeV) or H° and A (in the high mass region with ~ 120 GeV), as shown 
in Figs. 10 and 11 of Ref. 
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Table 2: Cross sections in fb for neutral Higgs boson production in the MSSM with tan/3 = 40, at the 
upgraded Tevatron and the LHC, are shown for four different CTEQ4 PDFs. They are separately 
given for the LO and NLO processes, and for the bb — > A°X and bg — ► A°X sub-processes. For the 
upgraded Tevatron the top number is for m\ = 200 GeV, the middle is for txia = 300 GeV, and the 
lowest is for ijia = 400 GeV. For the LHC the top number is for ijia = 400 GeV, the middle is for 
m,A = 700 GeV, and the lowest is for itla = 1 TeV. 
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Figure 12: LO and NLO cross sections for the neutral Higgs A production in the MSSM with 
tan/3 = 40, at the Tevatron and the LHC. (a) For each collider we show the NLO cross sections 
with the resummed running Yukawa coupling (solid) and with one- loop Yukawa coupling (dashed), 
as well as the LO cross sections with resummed running Yukawa coupling (dotted) and with tree- 
level Yukawa coupling (dash-dotted), (b) The NLO (solid), the bb (dashed) and bg (dash-dotted) 
sub-contributions, and the LO (dotted) contributions are shown. Since the bg cross sections are 
negative, they are multiplied by —1 in the plot. The cross sections at y/S = 1.8 TeV are multiplied 
by 0.1 to avoid overlap with the y/S = 2 TeV curves, (c) The NLO cross sections with QCD running 
Yukawa coupling (solid curves) and those with additional SUSY correction to the running coupling 
are shown (upper dashed lines for the Higgs-mixing parameter fi = +500 GeV and lower dashed 
lines for fi = -500 GeV). 
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Figure 13: The K-factors for the ^4° production in the MSSM with tan/? = 40 are shown for the 
NLO [K = <7nlo/o"lo, solid lines), bb (K = cr^/oxo = (clo + A<7 fc g) / f aho, dashed lines), and bg 
(K = — Ao"5 g /o"Lo ; dash-dotted lines) contributions, at the upgraded Tevatron (a) and the LHC 
(b). 



We then consider the large bottom Yukawa coupling of the neutral 6-Higgs (h®) and 6-pion 
(7r°) in the TopC model 0, [II], [33fl. The new strong U(l) force in this model is attractive 
in the <it> channel but repulsive in the <bb> channel. Thus, the top but not the bottom 
acquires dynamical mass from the vacuum. This makes the t- Yukawa coupling (y t ) super- 
critical while the b- Yukawa coupling (yb) sub-critical, at the TopC breaking scale A, i.e., 
j/b(A) ~ y cr it = a/87t 2 /3 ~ ?/t(A) , which requires yb being close to y t and thus naturally 
large. Our recent renormalization group analysis [R3| shows that the relation ?/fe(yu)~?/ i (/i) 
holds well at any scale fi below A. For the current numerical analysis, we shall choose a 
typical value of ?/ft(m t ) ~ ytimt) ~ 3, i.e., \C^\ = \C b ^\ ~ 3/y/2. In Fig. |TB|, we plot 
the production cross sections of h® or 7r° at the Tevatron and the LHC. This is similar to 
the charged top-pion production in Fig. f|, except the non-trivial differences in the Yukawa 
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Figure 14: The ratios of NLO cross sections computed by four different sets of CTEQ4 PDFs 
relative to that by the CTEQ4M for neutral ^-production in the MSSM with tan/3 = 40, at the 
upgraded Tevatron (a) and the LHC (b). 

couplings (due to the different tree-level values and the running behaviors) and the charm 
versus bottom parton luminosities. 



5. Conclusions 

In summary, we have presented the complete 0(a s ) QCD corrections to the charged scalar 
or pseudo-scalar production via the partonic heavy quark fusion process at hadron colliders. 
We found that the overall NLO corrections to the pp/pp — > ± processes are positive for 
above ~150 (200) GeV and lie below ~15 (10)% for the Tevatron (LHC) in the relevant range 
of (cf. Fig. H). The inclusion of the NLO contributions thus justifies and improves our 
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Figure 15: Transverse momentum distributions of pseudo-scalar A produced via hadronic colli- 
sions, calculated in the MSSM with tan/3 = 40. The resummed (solid) and 0(a s ) (dashed) curves 
are shown for uia = 200, 250, and 300 GeV at the upgraded Tevatron (a), and for uia = 250, 400, 
and 550 GeV at the LHC (b). 

recent LO analysis []!]. The uncertainties of the NLO rates due to the different PDFs are 
systematically examined and are found to be around 20% (cf. Table 1 and Fig. ^|). The QCD 
resummation to include the effects of multiple soft-gluon radiation is also performed, which 
provides a better prediction of the transverse momentum (Qt) distribution of the scalar 
<f) 0,± , and is important for extracting the experimental signals (cf. Fig. [?p. We find that the 
resummed total rate differs from the 0(a s ) rate only by a few percents which indicate that 
the size of the higher order corrections are likely much smaller than the uncertainty from 
the parton distribution functions (cf. Figs. 6 and 14). We confirm that the 2TeV Tevatron 
(with a 2 — 10 fb _1 integrated luminosity) is able to explore the natural mass range of the 
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m* (GeV) 

Figure 16: Cross sections for the neutral 6-pion ir® or 6-Higgs h® production via the bb- fusion in the 
TopC model at the Tevatron and the LHC. The NLO (solid), the qq' (dashed) and qg (dash-dotted) 
sub-contributions, and the LO (dotted) contributions with resummed running Yukawa coupling are 
shown. Since the qg cross sections are negative, they are multiplied by —1 in the plot. The cross 
sections at y/S =1.8 TeV are multiplied by 0.1 to avoid overlap with the V~S = 2 TeV curves. 

top-pions up to about 300-350 GeV in the TopC model 0, [IT]] for the typical tR-CR mixing 
of Ky R ~0.2 — 0.33 [cf. eq. ([11])]. Measuring the top polarization in the single-top event will 
further improve the signal identification. On the other hand, due to a possibly smaller (p^- 
b-c coupling in the 2HDM, we show that to probe the charged Higgs boson with mass above 
200 GeV in this model may require a high luminosity Tevatron (with a 10 — 30 fb _1 integrated 
luminosity). The LHC will further probe the charged Higgs boson of the 2HDM up to about 
0(1) TeV via the single-top and W^h (or W ± H°) production. The complementary roles of 
the tb and W ± h° channels in the different regions of the Higgs mass and the Higgs mixing 
angle a are demonstrated. We have also analyzed a direct extension of our NLO results to the 
neutral (pseudo-) scalar production via the 66-fusion for the neutral Higgs bosons (A , h°, H°) 
in the MSSM with large tan/3, and for the neutral 6-pion (tt®) or 6-Higgs (h®) in the TopC 
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model with U (l)-tilted large bottom Yukawa coupling. In comparison with the <p°hh associate 
production []33| , this inclusive 0°-production mechanism provides a complementary probe for 
a neutral Higgs boson (with relatively large mass), whose decay products, e.g., in the bb or rr 
channel, typically have high transverse momenta (~ and can be effectively detected 



35[ . This is particularly helpful for the discovery reach of the Tevatron. Further detailed 
Monte Carlo analyses at the detector level should be carried out to finally conclude the 
sensitivity of the Tevatron Run-II and the LHC via this process. 

At the final stage of writing up this manuscript, we became aware of a new preprint EC 



which studied the QCD corrections for the neutral Higgs production bb — > H° within the 
SM, and partially overlaps with our Sec. 4 as the pure NLO QCD correction is concerned. 
The overlapped part is in general agreement with ours except that we determine the counter 
term of the Yukawa coupling (expressed in terms of the relevant quark mass) by the on-shell 
scheme (cf. Refs. P7fl ) while Ref. [[40] used MS scheme. After resumming the leading 
logarithms into the running mass or Yukawa coupling, the two results coincide. Note that the 
apparent large 0(a s ) correction derived in Ref. [40] is due to the fact that it only includes the 



contribution from the bb sub-process, which is part of our complete 0(a s ) contribution. The 
inclusion of the NLO contribution from the bg sub-process, which turns out to be negative 
and partially cancels the bb contribution, yields a typical size of 0(a s ) correction to the 
production rate of a neutral Higgs boson produced via heavy quark fusion. The bg sub- 



process is identified as 0(1/ ln[m#/mft]) instead of 0(a s ) correction in Ref. |40 
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Appendix 



In this appendix, we present the individual NLO parton cross sections computed at 
D = 4 — 2e dimensions. We note that, unlike the usual Drell-Yan type processes, the one-loop 
virtual contributions (cf. Fig. lb-d) are not ultraviolet (UV) finite unless the new counter 
term from Yukawa coupling (related to the quark-mass renormalization, cf. Fig. le) is in- 
cluded. 



A. Partonic processes cb — ► <p + X 

The spin- and color-averaged amplitude-square for the cb 
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2nC f 



2 = (\Cl? 



\Cr\ 2 )^ 



(l-r)(i + | + 2 
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sm 2 
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The individual contributions (from the virtual loop and real gluon emission) to the NLO 
partonic cross section are: 
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where the standard plus prescription (•■■)+ is given by 



da X (a) [£(«)-£(!)]. 



(24) 



In (|23"D , the infrared — poles cancel between AcT^p Ual and Acx^ al . The term AcT[™p Ual from the 



3 3 
1 with e uv = -e IR 



virtual loop actually contains two types of - poles inside [• • 

3 

e = (4 — D)/2 > 0. Also, the — pole inside the Yukawa counter-term contribution Act 

3 e 

is ultraviolet while the H — pole inside AcT r | al is infrared (IR). We see that the contribution 

e 



virtual 
count 
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Account* 1 from the counter term of the Yukawa coupling is crucial for cancelling the UV 



divergence from A?^ (which is absent in the usual Drell-Yan type processes), while the 
soft - divergences between A?^™ 1 and Aa^ al cancel. Finally, the - collinear singularity 
inside A?™ 1 will be absorbed into the re-definition of the PDF via the quark-quark transition 
function Pq^-q(r). All the finite terms are summarized in Eq. (f|). 

B. Partonic processes gc,gb — > <p + X 

The spin- and color-averaged amplitude-square for the gc, gb — > <p + X process is 

,2 ' 



3(1 - e) 



(|c L | 2 + i^fV 6 




2 -2- 



u m A 



st 



(25) 



The 0(a s ) partonic cross section for the quark-gluon fusions is given by: 
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^(r) + 7(-3 + 7?)(l-f) 



(26) 



where it is clear that the collinear - singularity will be absorbed into the re-definition of 

e 

the PDF via the gluon-splitting function -PgH ff (r). The final result is finite and is given in 
Eq. (|). 
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